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We study the process e* e^ > п'я (3686) using 5.1 fb ^! of data collected at 16 center-of-mass 
energy (4/3) points from 4.008 to 4.600 GeV by the BESIII detector operating at the BEPCII collider. 
The measured Born cross sections for e* e^ — п'я (3686) are consistent with previous results, 
but with much improved precision. А fit to the cross section shows contributions from two structures: 
the first has M = 4209.5 + 7.4 + 1.4 MeV/c? and Г = 80.1 + 24.6 + 2.9 MeV, and the second has 
M = 4383.8+4.2+0.8 MeV/c? and T = 84.2+12.5+2.1 MeV, where the first errors are statistical and 
the second systematic. The lower-mass resonance is observed in the process ete” — п'я (3686) 
for the first time with a statistical significance of 5.87. A charged charmonium-like structure is 
observed in the т (3686) invariant mass spectrum for data аб Vs = 4.416 GeV. A fit with an S- 
wave Breit- Wigner function yields a mass M = 4032.1 3- 2.4 MeV/c’, where the errors are statistical 
only. However, there are still unresolved discrepancies between the fit model and data. The width 
of the intermediate state varies in a wide range for different kinematic regions within the data set. 
Therefore no simple interpretation of the data has been found, and a future data sample with larger 
statistics and more theoretical input will be required to better understand this issue. 


PACS numbers: 14.40.Rt, 14.40.Pq, 13.66.Bc 


І. INTRODUCTION structures. There are two resonances observed in the 
process етет — ntr” Ро in the mass region between 
Бе 4.2 and 4.4 GeV/c?. The lower mass state, У (4220), 

In the past decade, a series of new vector (/ © = аз а mass consistent with the lower mass state ob- 
177) charmoniumlike states, e.g., the Y (4260), Y (4360), served in ete” > ля Jib |J], but has а somewhat 
and Y (4660), have been observed in етет annihilation larger width. The higher mass state, У (4390), does not 
through dipion hadronic transitions to low mass char- match any known vector charmonium or charmonium- 
monium states [H5]. Many theoretical interpretations like states. Therefore, a more precise measurement of 
have been proposed to understand the underlying struc- the cross section of e+e- — mtr (3686) at BESIII 
ture of the Y family of states, such as hybrid charmoni- 


um [6], tetraguark |7], or hadronic molecule Е etc. The 

00 was first observed in Ше — ув У (4360) 一 In recent years, another mysterious and interesting new 
ysr T” TO (3686) by BaBar 4 and subsequently con- | pattern of charmoniumlike states, the Zz's, have been 
firmed by Belle à Recently, two resonant structures observed both in final states containing a «аква, ріоп 
were observed in the processes ete” — п'я J/w at and a low mass charmonium state |3 mi 113], and in pairs 
BESIII [9], indicating the i 4260) resonance reported оЁ charmed mesons (DC? D*)* |14116). А similar pattern 
by previous experiments actually consists of two of states is also observed in the bottomonium system ከሻ. 


will help to clarify the spectrum of vector particles. 
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All of these states contain a minimum of four quarks 
as they are electrically charged and their constituents 
include a cc pair. They are good candidates for meson- 
meson molecules or tightly bound tetraquark states. To 
clarify their true nature, it is essential to search for this 
kind of exotic state in other final states. Belle's recently 
updated result of етет > yrsr t” тт (3686), using their 
full data samples, shows evidence for а new 22 state 
in the mass spectrum of т. (3686) with a mass around 
4.05 GeV/c? and a statistical significance of 3.50 回 . It 
is important to confirm the existence of this new 22 
candidate in other experiments. 

In this paper, we present a study of e+e- 一 
тїл (3686), with two decay modes (3686) — 
теп J/w (mode I) and (3686) — neutrals + J/w 
(mode II), where ‘neutrals’ refers to поло, n°, p and 
yy, at center-of-mass (c.m.) energies Vs from 4.008 to 
4.600 GeV. The J/4 is reconstructed in its promi- 
nent decay mode J/y 一 £*£^ [ያ = e/u). The actu- 
al c.m. energies are measured by studying the process 
ете” — ши" [18]. The data samples used in this 
analysis were collected with the BESIII detector at the 
BEPCII collider at 16 different c.m. energy points with 
a total integrated luminosity of 5.1 fb! 


II. THE BESIII EXPERIMENT AND THE DATA 
SETS 


BEPCII is а double-ring e*e^ collider running at 
c.m. energies between 2.0 and 4.6 GeV, reaching a peak 
luminosity of 1.0 x 10??cm-?s-! at a c.m. energy of 
3770 MeV. The cylindrical BESIII detector has an ef- 
fective geometrical acceptance of 9396 of 4x and is divid- 
ed into a barrel section and two endcaps. It contains a 
small cell, helium-based (60% He, 40% Са На) main drift 
chamber (MDC) which provides a momentum measure- 
ment of charged particles with a resolution of 0.596 at 
а momentum of 1 GeV/c in a magnetic field of 1 Tesla. 
The energy loss measurement (dE/dx) provided by the 
MDC has a resolution better than 696. A time-of-flight 
system (TOF) consisting of 5-cm-thick plastic scintilla- 
tors can measure the flight time of charged particles with 
a time resolution of 80 ps in the barrel and 110 ps in the 
endcaps. An electromagnetic calorimeter (EMC) consist- 
ing of 6240 CsI (ТІ) in a cylindrical structure and two 
endcaps is used to measure the energies of photons and 
electrons. The energy resolution of the EMC is 2.596 in 
the barrel and 5.096 in the endcaps for photon/electrons 
with an energy of 1 GeV. The position resolution of the 
EMC is 6 mm in the barrel and 9 mm in the endcaps. A 
muon system (MUC) consisting of 1000 m? of Resistive 
Plate Chambers (RPC) is used to identify muons and 
provides a spatial resolution better than 2 cm. A de- 
tailed description of the BESIII detector can be found in 
Ref. 54. 

The GEANT4-based Monte Carlo (МС) simula- 
tion software package Boost [22] is used to generate 


the signal and background MC samples. The етет colli- 
sion is simulated with the KKMC generator taking 
into account the beam energy spread and Initial State 
Radiation (ISR), where the line shape of the e*e^ 一 
тя (3686) cross section is taken from the latest Belle 
results [5]. The QED Final State Radiative (FSR) correc- 
tion for J/ 一 £* £^ is incorporated with PHOTOS А 
Large signal MC samples етет —> тя (3686) as well 
as the dominant backgrounds are generated exclusive- 
ly at each c.m. energy, where the decay (3686) — 
Tt) т (9) J/ is simulated with the JPIPI model |25], 
and the decay ete“ — т?л” (3686) is generated with 
both the JPIPI model and a model distributed uniformly 
in phase space (PHSP). The JPIPI model is construct- 
ed with an effective Lagrangian for the hadronic decays 
and the coupling constants in the model are obtained 
by a fit to the invariant 1777 mass distributions in 
w (3686) 一 777" J/w. The generic inclusive MC samples 
at ላ/# = 4.258 and 4.358 GeV, with integrated luminosi- 
ty equivalent to that of data, are generated to study the 
potential background contributions, where the known de- 
cay modes are generated with EvTGEN [26] with branch- 
ing fractions set to values in the Particle Data Group 
(PDG) |27, and the remaining unknown ones are gener- 
ated with LUNDCHARM 1281. 


III. EVENT SELECTION 


The final states of the decay етет -» п'я (3686) 
have topologies п'я ፳7”፳ (^£. and тїл £6 + 
Му (М > 2) for mode I and II, respectively. А good 
charged track must have a polar angle | cos0| < 0.93 in 
the MDC, and have the point of closest approach to the 
interaction point within 1 cm in the plane perpendic- 
шаг to the beam and within +10 cm along the beam. 
Candidate events are required to have five (with +1 net 
charge) or six (with zero net charge) charged tracks for 
mode I and four (with zero net charge) charged tracks for 
mode II. Photon candidates are reconstructed by clus- 
tering EMC energy deposits. They must have a mini- 
mum energy of 25 MeV for | cos6| < 0.8 or 50 MeV for 
0.86 < |cos0| < 0.92, and be isolated from all charged 
tracks by an angle > 10°. EMC timing is used to sup- 
press electronic noise and energy deposition unrelated to 
the event. At least two photon candidates are required 
for mode II. Since pions and leptons are kinematically 
well separated for the signal process, charged tracks with 
momentum p > 1.0 GeV/c are assumed to be leptons, 
while those with p « 0.65 GeV/c are pions. Electron 
and muon separation is performed with the deposited 
energy E in the electromagnetic calorimeter. Electrons 
must have E/pc > 0.7, while muons have E < 0.45 GeV. 
Candidate events are required to have one lepton pair 
777”, and the other tracks must be pions. The number 
of hits in the muon counter is further required to be not 
less than 5 for J/v > ሀ” ሀ” candidates. 

In mode I, a four-constraint (4C) kinematic fit impos- 
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ing energy-momentum conservation under the етет 一 
atin ntn L E hypothesis is applied for the event can- 
didates with six charged tracks, and a one-constraint 
(1C) kinematic fit with one missing pion is applied for the 
candidates with five charged tracks. The corresponding 
x? of the fit must be х2с < 60 or Xic < 15, respectively. 
The J/w signal appears with low background in the spec- 
trum of the ¢*é~ invariant mass, M(£* £^), calculated 
with the corrected momentum after the kinematic fit, and 
is required to be within 3.05 < M(£*£^) < 3.15 GeV/c?. 
To improve the mass resolution of the /(3686), the select- 
ed candidates with six or five good charged tracks are fur- 
ther fed into a five-constraint (5C) or two-constraint (2C) 
kinematic fit with an additional constraint on M (£*£^) 
to the nominal J/w mass (M(J/w)) 23, respectively. 
The (3686) signal is reconstructed with the 1777 J/w 
system (four combinations per event) whose invariant 
mass М(п'т J/1)), calculated with the corrected mo- 
mentum after the 5C or 2C fit, is closest to the nominal 
%(3686) mass (M (v(3686))) [2]. 

In mode II, no kinematic fit is applied. The J/wv 
signal is extracted by requiring 3.05 < M(t) < 
3.15 GeV/c?. Тһе (3686) signal is reconstruct- 
ed using the mass recoiling against the TTT“ sys- 
tem, Moil (atr). The requirement |М<%%6) 一 
M(w(3686))| > 8 MeV/c? is imposed to veto the 
background e*e^ — neutrals + (3686) (((3686) 一 
тїт Jh) "where Мо%ө) = М(тїт ያያን 一 
M(£*07) + MŲJ/4). A requirement |[ዳ/(ግግጽ?፳ ) – 
M(n)| > 50 MeV/c? is used to eliminate the background 
ете” => nJ/vV (n — тіт 79), where уу are the two 
photons with largest energy, and M(n) is the nominal 
7] mass 5ሻ. А requirement сов0,-л- < 0.9 is also ap- 
plied to remove the radiative Bhabha апа dimuon back- 
ground in which a photon converts into a ር” ር” pair and 
is misidentified as a 7777 pair. 


IV. EXTRACTION OF THE BORN CROSS 
SECTION 


With the above selection criteria, a prominent (3686) 
signal over a small background is observed in the 
М(п п" //ህ) and М" (л+л) spectra for Mode I and 
Mode II, respectively, as shown in Fig. [| for data at 
8 = 4.416 GeV as an example. Potential backgrounds 
without a J/w in the final state are explored with events 
within the J/w sideband regions. Backgrounds including 
a J/w in the final state are studied with both inclusive 
and exclusive MC samples. From these studies, no peak- 
ing background is expected in the range of interest in 
the two (3686) decay modes. To determine the signal 
yields, unbinned maximum likelihood fits are performed 
to the M(nx*«- 7/ф) and М" (л+л) spectra for the 
two modes, respectively. Here, the signal probability den- 
sity function (PDF) is described with the MC simulated 
shape convoluted with a Gaussian function to account 
for the mass resolution difference between data and MC 


simulation, and the PDF for background is presented as 
a linear function. Simultaneous fits to data at 16 differ- 
ent c.m. energies are performed. The c.m. energies share 
a common background shape and a common data/MC 
resolution difference. 
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FIG. 1. Distributions of (left) M(r+r J/w) for mode I and 
(right) እያ”) for mode II at Vs =4.416 GeV, with 
the integrated luminosity of 1074 pb !. Dots with error bars 
are data, the shaded histograms are background from J/w 
sideband, and the curves are the fit described in the text. 


The Born cross section is determined from 


А Nobs 
(р, 9 1 
i fur (+ô) (1109) B5 (i Й ), ( ) 


where eee is the signal yield extracted from the fit, Lint 
is the integrated luminosity, 1 + 0" is the ISR correction 
factor, 1 + ô” is the vacuum polarization factor, Br is 
the product of the decay branching fractions of the in- 
termediate states in the cascade decay (4.11% for mode 
I and 2.95% for mode II), and є is the detection efficien- 
cy. The various numbers used in the Born cross section 
calculation are summarized in Table [I] 

Several sources of systematic uncertainty are consid- 
ered in the measurement of the Born cross section. The 
uncertainty in the tracking efficiency for leptons and pio- 
ns are 1.0% per track 50. The uncertainty in the photon 
efficiency is 1.0% per photon [31]. The uncertainty relat- 
ed to the reguirement on the number of hits in the muon 
counter is 4.2%, determined by studying a control sample 
of етет + пп Ј/р. The uncertainty associated with 
the kinematic fit is estimated with the same technique 
as in Ref. [23]. The uncertainty associated with the J/4 
mass requirement is estimated by smearing the M (£*£^) 
distribution of MC samples according to the resolution 
of M(£*£^) in data, and the difference in MC efficien- 
cies is taken as the uncertainty. In mode II, the uncer- 
tainties associated with the requirements on М 0(3686)› 
М(үүт*т) and cos0,-,- are evaluated by changing 
the requirements. The uncertainty associated with the 
fit procedure is investigated by varying the fit range, re- 
placing the linear function by a second-order polynomial 
function for the background, and varying the width of the 
Gaussian function for signal. For data sets with large lu- 
minosity, the detection efficiency is estimated with a MC 
sample weighted according to the observed Dalitz plot 
М?(770(3686)) versus M?(x*--). The corresponding 
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TABLE I. Summary of the measurement of the Born cross section c? at individual c. m. energies. The subscript 1 or 2 
denotes mode I or II, respectively. The first uncertainties are statistical, and the second systematic. An upper limit at the 90% 
confidence level (C.L.) is determined by a profile likelihood method [29] for data samples with low signal significance. 


уз (GeV) 
- 0. Е 2. 4. ; ; 


4.008 0.0 + 0.6 0.2 3 
4.085 4.0 + 2.0 了 
4.189 3.8 + 2.0 引 
4.208 8.9 + 3.3 了 
4.217 13.0 + 3.6 £ 
4.226 315 + 18 

4.242 11.0 3- 5.9 

4.258 241 + 16 

4.308 : 17.0+ 4.2 

4.358 439 + 21 

4.387 56.6 + 7.6 

4.416 693 4 

4.467 15.1 3- 4.2 

4.527 13.4 + 4.0 

4.575 45 d- 2.3 

4.600 106 4 


uncertainty is estimated by varying the weight factors ac- 
cording to the statistical uncertainty of data in each bin 
on the Dalitz plane. For data sets with low luminosity, 
the detection efficiency is estimated using MC samples 
with the JPIPI model; the corresponding uncertainty is 
evaluated with large luminosity data sets, by taking the 
largest difference of efficiencies between the JPIPI mod- 
el MC sample and their nominal values with weighted 
MC. The uncertainty related with the ISR correction fac- 
tor is studied by replacing the input cross section line 
shape with the latest results from ВаВаг [1 in ККМС, 
and the change оп (1 + 6”) є is taken as an uncertain- 
ty. The uncertainty in the vacuum polarization factor 
is 0.5% taken from a QED calculation [33]. The uncer- 
tainty on the integrated luminosity is 1%, determined 
with large angle Bhabha events The uncertainties 
in the decay branching fractions of intermediate states 
are quoted from the PDG 54. The uncertainty from 
others sources, such as lepton separation, trigger efficien- 
cy, and FSR are negligible, and are conservatively taken 
to be 1.0%. Assuming all sources of systematic are inde- 
pendent, the total uncertainties are obtained by adding 
the individual values in quadrature, and are in a range 
between 7.7% to 14.1% and 7.4% to 20.1%, depending on 
c.m. energy, for mode I and II, respectively. 


The measured Born cross sections of e*e^ 一 
тїт (3686) at individual c.m. energies for the two 
(3686) decay modes are consistent with each other with- 
in their uncertainties. The measurements are there- 
fore combined by considering the correlated and un- 
correlated uncertainties between the two modes, accord- 
ing to Refs. |34, [35]. The comparison of {һе combined 
Born cross section of e* e^ — тт (3686) with those 
from previous experimental results is shown in Fig. [2] 
The results are consistent with former experiments, and 
have much improved precision. 


0.70 
0.75 


0.76 


Æ BESIII 
—- Belle 


ak 
о 
о 


æ BaBar 


Cross Section (pb) 
5 መ መ 


го 
о 


о 


710 42 43 44 45 46 
(s(GeV) 


FIG. 2. Born cross section of e*e^ 一 тїт w(3686). The 
dots (red) are the results obtained in this analysis, the trian- 
gles (green) and squares (blue) are from BELLE and Babar's 
latest updated results, respectively. The solid curve is the fit 
to BESIII results with the coherent sum of three Breit- Wigner 
functions. The dashed curve (pink) is the fit to BESIII results 
with the coherent sum of two Breit- Wigner functions without 
the Y (4220) hypothesis. The arrows mark the locations of 
four energy points with large luminosities. 


V. FIT TO THE CROSS SECTION 


To study the possible resonant structures in e* e^ 一 
тїт (3686), a binned x? fit is applied to describe the 
cross section obtained in this analysis in a energy region 
from 4.085 GeV to 4.600 GeV. Assuming that three res- 
onances exist, the PDF can be parameterized as 


A= ће + fa + fae, (2) 


where f; is for the Y (4220), / is for the Y (4390), fs is 
for the Y (4660), фі is the phase angle between Y (4390) 
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and Y (4220), and фә is the phase angle between У (4390) 
and Y (4660). 

The amplitude f; (i = 1,2,3) for each resonance is a 
P-wave Breit-Wigner function, defined as 


M; V1l2rBiT Г; faM) 
f= м мі- M2 +iM,T; 9(M;)' 3) 


where M; is the mass of the resonance, pee is the par- 
tial width бо є" e^, Г, is the total width, B; is the branch- 
ing fraction of the resonance decaying to T" 7" (3686), 
and Ф(М;) is the three-body PHSP factor. In the fit, M;, 
I;,the product Bee, and the relative phase angle 6; 
(1 = 1,2) are free parameters. The parameters of the 
Y (4660) are fixed to Belle's latest results 回 . 

The fit result is shown in Fig. Blas the solid curve and 
the resulting parameters of the resonances are summa- 
rized in Table There are two solutions with equally 
good fit quality, X?/ndf = 9.98/7, where паї is the num- 
ber of degrees of freedom. An alternative fit without 
Y (4220) is performed as shown by the dashed curve in 
Fig. D] yielding x2/ndf = 54.54/11. The significance of 
the Y (4220) is determined to be 5.80. 


TABLE II. Results of the fit to the ete” — пт (3686) 
cross section. The error is statistical only. 


Parameters Solution I Solution II 

M(Y4220) (MeV/c?) 4209.5 ፲ 7.4 
Г(У(4220)) (Мем) 80.1 + 24.6 

BTe * (Y(4220) (eV) 0.8 + 0.7 0.4 3 0.3 
M(Y4390) (MeV/c?) 4383.8 + 4.2 
Г(У(4390)) (Мем) 84.2 + 12.5 

Вге'є (Y(4390) (eV) 3.6 + 1.5 2.7+1.0 

фі (rad) 3.3 ፳ 1.0 2.8 + 0.4 

фә (rad) 0.8 + 0.9 4.7+0.1 


The systematic uncertainties accociated with the reso- 
nance parameters include the following. The uncertainty 
from the c.m. energy measurement is studied by takin 
the uncertainty of the c.m. energy A(./s)=0.8 GeV 
into consideration. The uncertainty due to the c.m. en- 
ergy spread is estimated by convoluting the fit PDF with 
a Gaussian function with a width of 1.6 MeV, which 
is the energy spread determined by the Beam Energy 
Measurement System [36]. The uncertainty from the 
fit range is investigated by including the first energy 
point Vs = 4.008 GeV or excluding the energy point 
Vs = 4.085 GeV. The uncertainty from the У (4660) 
resonance is studied by varying the parameters of the 
У (4660) within its uncertainties. Table [ITI] summarizes 
the systematic uncertainties for the mass and total width 
of the resonances Y (4220) and Y (4390). 

A comparison of the masses and widths of the Y (4220) 
and Y (4390) determined in this work to those of other 
vector states observed at BESIII in the mass region be- 
tween 4.2 and 4.4 GeV/c? are shown in Table [IV] 


TABLE III. Systematic uncertainty in the measurement of 
resonance parameters. 


Sources Y (4220) Y (4390) 
M (MeV/c?) T (MeV) M (MeV/c?) T (MeV) 
Energy scale 0.1 0.0 0.1 0.0 
Energy spread 1.3 0.1 0.7 1.6 
Fit range 0.2 1.8 0.2 0.5 
У (4660) 0.6 2.3 0.3 1.2 
"Total 14 2.9 0.8 2.1 


VI. STUDY OF INTERMEDIATE STATES 


Intermediate states in the process ete  — 
тт (3686) are investigated in the data samples 
that have large integrated luminosity. A requirement of 
3.68 < М(тіл”.УФ)(Мтезі(д%д-)) < 3.70 GeV/c? 
is applied to extract the (3686) signal, and the side- 
band regions, 3.63 < М(п' п J/A)(Mrecll(g t 67)) < 
3.65 GeV/c? or 3.73 < М(л+л- //ф) < 3.75 GeV/c?, 
are used to explore the potential поп-/(3686) back- 
grounds, where only the left side sideband region is 
used in mode II since a long tail appears on the right 
side of the (3686) signal due to ISR. The поп- (3686) 
backgrounds are found to be small, and do not produce 
peaks in the various distributions. 

With the above selection criteria, the Dalitz plots of 
Мп п" ) versus Mž(1=4(3686)) and the correspond- 
ing one-dimensional projections are shown in Fig. [8] for 
data samples at 1/8 = 4.416, 4.358, 4.258 and 4.226 GeV, 
individually, where the plots include the candidates of 
the two (3686) decay modes. It should be noted that 
the two higher energy points are from the Y (4390) peak 
region, while the other two are from the Y (4220) peak re- 
gion. These four points are marked with vertical arrows 
in Fig. 

For data at ys = 4.416 GeV, a prominent nar- 
row structure is observed around 4030 MeV/c? in the 
M (1—15(3686)) spectrum. The structure is also evident 
in the corresponding Dalitz plot, but it appears to be 
more complex when looking at high and low M(z*7-) 
ranges. In the low М(т"т ) region, there seem to be 
two separate т" (3686) structures, presumably corre- 
sponding to a physical structure and its kinematic re- 
flection. But in the high М(т" т ) region, only one 
broad т" (3686) structure is seen. For data at Vs = 
4.358 GeV, there is no obvious structure observed in 
the M (1 —715(3686)) spectrum, but a cluster of events ap- 
pear in the low М(т'"т ) region on the corresponding 
Dalitz plot. It is worth noting that, at this c.m. en- 
ergy, a physical structure with a mass of 4030 MeV/c? 
in the M(1*4(3686)) spectrum has a reflection at the 
same mass position. For data at Vs = 4.258 GeV, two 
bumps around 3900 and 4030 MeV/c? are visible in both 
the Dalitz plot and in the M (x-^45(3686)) spectrum. For 
data at ys = 4.258 GeV, the possible structures with 
masses of 3900 and 4030 MeV/c? in the M (7>4(3686)) 


TABLE ТУ. Masses апа widths ої the vector charmonium states observed from different processes at BESIII in the mass region 
between 4.2 and 4.4 GeV/c”. The subscript 1 ог 2 denotes the lower mass state or higher mass state. 


Process 
ете” > ЧХсо 
ete” a ntr J/y 
ete- 一 ntr he 
ete- - nt DU D*- + с.с 
ете” — mm (3686) 


Му (MeV/c) T; (MeV) М» (MeV/c?) Гә (MeV) 


10.9 7.0 | 101.4*25 + 10.2[9] 
十 5.5 12.3 十 6.3 十 16.2 
4218.4173 4:09. 66.07, E04 43591579310 — 139.5+203 士 0.6 [10] 
4224.8 +5.6+4.0 723491409 4400.12 9.3: 2.1 — 181.7 + 16.9 + 7.4 [38] 
4209.5 d- 7.5 1.4 | 80.1- 24.63: 2.9 | 4383.8 + 4.2 + 0.8 84.2 + 12.5 + 2.1 
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FIG. 3. Dalitz plots of М?(т<)(3686)) versus M? (rtr), distributions of M?(z^1»(3686)) (two entries per 


event), and M?(x*«-) for data at /5 = 4.416, 4.358, 4.258 and 4.226 GeV, with integrated luminosities of 1074, 
540, 826 and 1092 pb“! respectively. Dots with errors are data. For the plots at \/s = 4.416, 4.358 and 4.258 GeV, 
the solid curves (red) are projections from the fit; the dashed curves (pink) show the shape of the intermediate 
state; the dash-dotted curves (blue) show the shape from the direct process ете” — т'т (3686) obtained 
from the JPIPI MC model; the shaded histograms (green) show the non-1/(3686) background estimated with the 
(3686) sideband. For plots at Vs = 4.226 GeV, the dashed (pink) and dash-dotted (blue) curves show the shapes 
from the intermediate state and the direct process ете — тіл (3686) (with arbitrary scale). In all plots, the 


two (3686) decay modes are combined. 
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spectrum have kinematic reflections at each other's mass 
positions. For data at \/s = 4.226 GeV, no structure is 
clearly seen, which is very different from the behavior at 
the energy point close by, ч/з = 4.258 GeV. A further 
striking feature for data at Vs = 4.226 GeV is a very 
different M (r^) distribution from those at the other 
three energy points. 

To characterize the structure observed on the 
M (x*~(3686)) spectrum for data at ላ/# = 4.416 GeV, 
an unbinned maximum likelihood fit is carried out on the 
Dalitz plot of M?(z*4(3686)) versus M?(x (3686)) 
(denoted as х and y in formula 4). Assuming an in- 
termediate state with spin parity 1*, the Dalitz plot is 
parameterized by the coherent sum of the process with 
an intermediate state and the direct process e*e^ — 
ля (3686). The PDF of the intermediate state is 
described with an S-wave Breit-Wigner function with- 
out considering interference among the charged conjugate 
modes, 


»-9/ሯ 
(M2 — y + М2 -T?/c*’ 
(4) 
where р (q) is the (3686) (intermediate state) momen- 
tum in the т (3686) (initial етет) rest frame, and Мр 
and Г are the mass and width of the intermediate state. 
The 2-dimensional mass resolution and the detection effi- 
ciency, determined from MC simulations, are incorporat- 
ed in the PDF for the intermediate states in the fit. The 
PDF ої the direct process етет — пт (3686) is taken 
from a MC-simulated shape using the JPIPI model, and 
that of the non- (3686) background is described with the 
distribution of events in the (3686) sideband region. A 
simultaneous fit constraining the mass and width of the 
intermediate state is carried out by maximize the product 
of the likelihood values of the two (3686) decay modes. 
The fit process is validated using MC samples. The dis- 
tributions of mass resolution and detection efficiency are 
provided in the appendix. 

The fit yields a mass of M = 4032.1 + 2.4 MeV/c? 
and a width of T = 26.1 + 5.3 MeV for the intermediate 
state with a significance of 9.20, evaluated by comparing 
the likelihood values with or without the intermediate 
states included. The fit projections оп M?(z+(3686)) 
and 37*(፳? п ) for data at Vs = 4.416 GeV are shown 
in Fig. [8] It can be seen that the overall fit curve does 
not match the peaking structure on the M (x^w(3686)) 
spectrum in data, and the corresponding confidence lev- 
el (C.L.) of the fit is only 896, as estimated by toy-MC 
tests. Alternative fits with different assumptions of the 
spin parity of the intermediate state, including the in- 
terference among the charge conjugated modes, and in- 
cluding the contribution of 2.(3900)= are explored. In 
these fits the parameters of the intermediate state are 
close to the norminal fit result and the fit qualities are 
not improved significantly. As shown in the Dalitz plot, 
the behavior of the structure is very different between 
the high M(z^7- ) region and the low M (^7) region. 


рес 


(M2 — a)? МА Te 


A similar fit to data with the additional requirement 
M? (ntr) > 0.3 (GeV/c?)? is performed, which yields 
a mass of M = 4030.3 + 0.1 MeV/c? and a width of 
Г = 5.1 +0.2 MeV. The corresponding projection of the 
fit and data on the М?(т:-7(3686)) distribution is shown 
in Fig. Д and the fit C.L. is 50%. 
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FIG. 4. Projection of M?(x-4(3686)) at VS = 4.416 GeV 
with a requirement of М?2(л"т”) > 0.3 (GeV/c?)?. 


Similar fits are carried out to data at ys = 4.358 
and 4.258 GeV, where the parameters of the interme- 
diate state are fixed to those from the fit to data at 
Vs = 4416 GeV. The projections оп М?2(7(3686)) 
and M?(x*«-) are also shown in Fig. B] The statistical 
significance of the intermediate state is 3.60 and 9.60 for 
data at ys = 4.358 and 4.258 GeV, respectively. For 
data at Vs = 4.358 GeV, as shown in the Dalitz plot, 
a cluster of events appears іп the М?(л:-)(3686)) spec- 
trum at low M?(z*«-), which also shows a hint of an 
intermediate structure. As mentioned previously, for da- 
ta at Vs — 4.258 GeV, the structures with masses of 
3900 and 4030 MeV/c? on the М(л:-0(3686)) spectrum 
are kinematic reflections of each other, so the fit results 
are strongly dependent on whether the Z.(3900)- is in- 
cluded in the fit or not. For data at Vs — 4.226 GeV, 
an intermediate state with a mass of 4030 MeV and its 
reflection are very close to the kinematic boundary of 
the three-body decay, so no obvious peak is observed in 
the M(z-—v(3686)) spectrum. The anomalous distribu- 
tion on the M (z^ п ) spectrum is, however, discussed in 
Ref. [39]. For the other energy points with high statis- 
tics, such as \/s = 4.387 and 4.600 GeV, the Dalitz plots 
and the distribution of M?(r*4(3686)) and M? (ntr) 
are shown in Fig. [5] 


VII. SUMMARY 


In summary, based on 5.1 Ұр”! of e+e- collision data 
with c.m. energies between 4.008 and 4.600 GeV, col- 
lected with the BESIII detector, the Born cross sec- 
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Dalitz plots of M2(rtwW%(3686)) versus M° (rtr), distributions of M?(m*7)(3686)) (two entries per event), and 


M*(n*1") for data аб 4/5 = 4.38 and 4.600 GeV, Dots with errors are data, the dashed (pink) and dash-dotted (blue) curves 
show the shapes from intermediate state and direct process e* e^ — т?л” w(3686) (with arbitrary scale). 


tions of етет — п'я (3686) are measured. The mea- 
sured cross sections are in good agreement with those 
from Babar and Belle [4 8, but with significantly im- 
proved precision. The cross section is fit with three co- 
herent Breit-Wigner functions. The parameters of two 
resonances, the Y(4220) апа Y(4390), are determined 
from the fit to be M = 4209.5 + 7.4 + 1.4 MeV/c? and 
T = 80.1 + 24.6 + 2.9 MeV for the Y (4220), and M = 
4383.8 + 4.2+0.8 MeV/c? and T = 84.2 2.1 MeV 


፡ 12.5 
for the Y(4390), where the first errors are statistical 
and the second systematic. The resonance У (4220) is 
observed in the process e*e^ — п'я (3686) for the 
first time with a significance of 5.80. Both of the reso- 
nances Y (4220) and Y (4390) are consistent with the vec- 
tor resonances observed іп ете” — ля አሩ [19]. The 
Y (4220) state is possibly the same state as observed in 
e+e- => ntr JI ала ete~ > охоо |37, and an at- 
tempt was made with а combined fit to the ር” e^ — шхо, 
ata Мф, ፳?፳” he, and я" DŪD*- + c.c cross sections 
measured by ВЕ5Ш to determine the resonant parame- 
ters of the У (4220) [49]. Peaking at а similar mass was 
also observed in the process етет — nJ/w at BESIII. 


A charged charmonium-like structure is observed in the 
M (п Ф(3686)) spectrum for data at Vs = 4.416 GeV. 
A fit to the structure, assuming the existence of a spin- 
parity 1* charmonium-like state, yields a mass M — 
4032.14+2.4 MeV /с?, where the errors are statistical only. 
However, the fit cannot describe data well. The width of 
the intermediate state varies over a wide range for dif- 
ferent kinematic regions of data at ys = 4.416 GeV. 
Similar fits are carried out to data at Vs — 4.258 and 
4.358 GeV, where there is also evidence for the new struc- 
ture at 4030 MeV/c? in the M(1*4(3686)) spectra. For 
data at ys = 4.226 GeV, no obvious structure in the 
M(x*w(3686)) spectrum is observed, but the M (ntr) 


spectrum shows an anomalous distribution. The mea- 
sured mass of the intermediate state deviates from that of 
the structure observed by Belle with their full data sam- 
ple [5] by over 30. Larger data sets, as well as additional 
theoretical input, are necessary for further understanding 
of the intermediate structures in ር? ይ” — пт v(3686). 
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FIG. 6. Resolution function for M?(x^45(3686)) at (a) 4.416, (b) 4.358, (c) 4.258 and (d) 4.226 GeV. The upper 
plots for mode I and bottom for mode II. Dots with error bars are data, the red solid line are fit results, dot-dashed 
line in pink denotes contribution from Gaussian function one, and the dashed line in blue denotes contribution 


from Gaussian function two. 
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